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1.  INTRODUCTION 


The  development  and  the  design  of  advanced  hypersonic  vehicles  such  as  the  Na¬ 
tional  Aerospace  Plane  (NASP)  or  the  Space  Shuttle  will  strongly  depend  on  computer 
simulation  mainly  because  some  of  the  flight  conditions  experienced  by  these  vehicles  will 
be  difficult  or  prohibitively  expensive  to  create  in  ground  facilities.  Computer  simulation 
will  be  used  to  reduce  the  need  for  expensive  test  time  or  will  replace  those  experimental 
tests  which  can  not  be  performed.  However,  computer  simulations  require  validation 
before  their  results  can  be  accepted.  These  validation  tests,  which  must  be  performed  in 
hypersonic  and  turbulent  air  flows,  require  the  simultaneous  measurement  of  air  tempera¬ 
ture,  density  and  velocity  at  a  high  temporal  and  spatial  resolution  which  will  allow  to 
resolve  turbulence  structures.  The  accuracy  of  the  diagnostic  system  should  be  in  the  neigh- 

\,  .  r  r-  ’  ■  ■ 

borhood  of  1%  and  the  system  must  respond  well  both  to  cold-air  flows  (>100  K)  in  the 
free  stream  and  hot-air  flows  (>1000  K)  experienced  at  the  stagnation  point.  Furthermore, 
elimination  of  air  seeding  requirements  will  simplify  the  measurement  and  will  allow  its 
application  in  any  facility. 

We  have  been  engaged  for  more  than  two  year  (under  NASA  contract  #  NAS2- 
12410)  in  the  development  of  a  technique  for  the  measurement  of  air  temperaiure  and  den¬ 
sity  using  an  Argon-Fluoride  (ArF)  laser.  The  technique  applies  Laser-Induced  Fluores¬ 
cence  (LIF)  in  molecular  oxygen  for  the  temperature  measurement,  and  Raman  scattering 
from  nitrogen  and  oxygen  for  the  density  measurement.  Our  efforts  have  already  identified 
the  measurement  potential  and  restrictions  [1,2]  of  the  proposed  technique  and  are  now 
directed  toward  the  demonstration  of  this  approach  and  towards  the  engineering  of  the 
working  device. 
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In  addition,  we  have  started  under  an  SBIR  contract  from  DARPA  [3]  and  under 
the  present  funding  to  investigate  the  possibility  of  applying  an  ArF  laser  for  velocity 
measurements  of  unseeded,  hypersonic,  cold-air  flows.  We  expect  that  this  extension  of  our 
work  will  lead  to  the  development  of  an  integrated  system  which  will  be  capable  of  measur¬ 
ing  simultaneously  the  temperature,  density  and  velocity  of  unseeded-air  flows. 

We  proposed  to  study  under  the  present  contract  two  velocity  measurement  tech¬ 
niques  which  will  both  rely  for  the  measurement  on  an  ArF  laser.  In  the  first  approach  the 
air  velocity  is  to  be  measured  by  using  an  intense  ArF  laser  pulse  which  dissociates  atmos¬ 
pheric  C02  into  CO  and  atomic  carbon.  The  atomic  carbon  is  a  strong  absorber  at  193. lnm 
and  a  strong  radiator  but  it  is  short  lived  radical  in  most  flow  situations.  The  CO  molecule 
on  the  other  hand  is  chemically  stable  and  may  reside  indefinitely  in  the  air  flow.  There¬ 
fore,  we  proposed  to  use  the  CO  molecule  as  a  tracer  in  a  modified  ’Time-of-Flight" 
velocimetry  technique.  We  have  already  identified  the  optimal  measurement  approach  [3], 
Here  we  investigated  the  effects  of  air  temperature  and  density  and  the  irregularities  in  the 
distribution  of  the  energy  across  the  laser  beam  on  the  accuracy  of  the  velocity  measure¬ 
ment.  We  have  also  conducted  a  preliminary  experiment  where  the  sensitivity  of  the  tech¬ 
nique  to  supersonic  air  flows  was  tested.  The  results  of  these  studies  are  reported  in  section 
2. 

In  the  second  approach  the  air  velocity  is  to  be  measured  by  the  Doppler  shift  ex¬ 
perienced  by  the  moving  absorbing  molecules.  In  this  technique  the  spectral  shift  between 
the  absorption  wavelength  of  molecules  at  rest  and  the  absorption  wavelength  of  the 
moving  molecule  is  measured.  We  have  already  estimated  the  potential  and  accuracy  of 
this  technique  [3].  Here,  we  evaluated  the  potential  of  this  technique  and  estimated  the  ef¬ 
fects  of  air  temperature  and  density  on  the  measurement  accuracy.  Since  the  signal  varies 
linearly  with  laser  intensity  we  did  not  have  to  investigate  the  effect  of  the  energy  distribu¬ 
tion  across  the  beam.  We  described  these  results  in  section  3. 
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2.  FLOW  TAGGING. WmJ  .CAJRB.OK  ATOMS 


The  C02  molecules  which  are  present  in  air  (at  .03%  mole  fraction)  can  be  dis¬ 
sociated  by  the  ArF  laser  into  CO  and  O  (see  figure  1).  At  193  nm  this  is  a  single-photon 
spin-forbidden  dissociation  into  ground-state  products  [4],  The  variation  of  the  absorption- 
coefficient  with  wavelength  in  the  UV  was  measured  [4]  in  a  cell  containing  700  torr  of  C02 
at  300  K.  It  was  shown  to  be  a  monotonously  decreasing  curve  with  an  absorption  coeffi¬ 
cient  at  193  nm  of  (^=0.1  m*1.  The  line-width  of  the  ArF  laser  is  narrow  relative  to  the 
spectral  variations  of  this  curve  and  therefore  this  absorption  coefficient  is  approximately 
constant  throughout  the  tuning  range  of  the  laser.  Although  this  is  a  weak  absorption,  a 
large  fraction  of  molecules  may  be  dissociated  when  irradiated  by  an  intense  ArF  laser. 
The  number  density  of  C02  molecules  at  700  torr  and  300  K  is  =  2.2T025  #/m3.  The 
absorption  cross-section  is  therefore: 

ai/Nco2  =  0.1/2, 2T025  =  4.55T0'27m2  (1) 

The  product  molecules  are  further  dissociated,  via  a  two-  photon  process,  into 
atomic  carbon  which  is  left  in  the  2^  state  [5].  The  two-photon  absorption  coefficient  for 
this  process  is  [6]  q2  =  5  TO* 11  m4/watt.  Finally,  the  carbon  atom  is  resonant  with  the  laser 
radiation  at  193.1  am  (see  figure  1)  and  may  therefore  be  excited  into  the  31?  state  from 
which  it  radiates  at  247.8  nm  while  decaying  into  the  21S  state  [5].  The  Einstein  A  coeffi¬ 
cient  for  the  transition  from  3^  to  21D  is  [7]  A^  =  2.4T08  sec'1  and  for  the  transition  from 
31?  to  2’S  is  [7]  ApS  =  3.4TQ7  sec'1.  While  the  interaction  of  CO  with  the  ArF  laser  was 
observed  by  Bokor,  et  al.[5]  the  overall  process  which  starts  with  atmospheric  COz  was  ob¬ 
served  recently  by  the  present  principal  investigator  [1]. 
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MULTI-PHOTON  INTERACTION  OF  ArF  LASER 

WITH  CO2 


C02  +  hv - ►CO  +  0 

CO  +  2hv - ►C(21D)  +  Q 

C(21D)  +hv-^C(31P) 

COV) - ►  C(21P)  +  M267.8  nm) 


Figure  1:  The  photo-dissociation  of  C02  and  CO  and  the  energy  level  diagram  for  the 

LIF  in  atom  carbon. 
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In  a  previous  work  [3]  we  determined  which  of  the  intermediate  products  of  this 
process  is  long-lived  and  therefore  the  most  likely  to  follow  the  flow  and  be  detected  at  a 
latter  time  by  delayed  probing.  It  was  demonstrated  that  carbon  atoms  in  the  C^D)  state 
are  short-lived.  They  recombine  with  N2  and  Oz  to  form  new  molecular  species.  Therefore, 
the  time  dependent  concentration  of  these  laser  generated  atoms  follows  very  closely  the 
temporal  variations  of  the  laser  pulse  itself.  Unless  the  flow  is  extremely  fast,  these  atoms 
will  combine  with  other  air  constituents  before  moving  out  of  the  laser  beam.  The  number 
density  of  the  CO  molecules  was  shown  to  remain  constant  long  after  the  laser  intensity  has 
diminished.  Thus,  only  the  CO  molecules,  which  are  stable,  may  be  generated  by  the  laser 
pulse  at  one  location  and  be  probed  downstream  to  determine  their  new  location  and 
hence  their  speed.  In  contrast,  the  carbon  atoms  which  are  eliminated  from  the  flow  by 
several  high-rate  mechanisms  may  not  serve  as  a  flow-tagging  species. 

2.1  SINGLE  PULSE  CARBON  VELOCIMETRY 

Although  the  generation  of  the  CO  molecule  by  the  ArF  laser  is  a  straightforward 
process,  the  probing  of  this  molecule  is  complex.  The  wavelength  of  the  dipole-allowed 
transition  between  the  ground-state  and  the  nearest  electronic  state  in  CO  is  in  the  vacuum 
UV  [8].  Thus,  single-photon  excitation  of  the  molecule  cannot  be  obtained  with  conven¬ 
tional  sources.  Instead,  a  two-photon  pumping  scheme  must  be  employed.  One  such 
scheme  was  already  demonstrated  for  trace  species  diagnostics  in  flames  [9].  In  this  tech¬ 
nique,  two  photons  at  230.1  nm  are  used  to  excite  the  CO  molecule  from  the  X1S+(v"=0) 
to  the  B1e+(v’=0)  state.  The  resulting  fluorescence  is  in  the  ETiT  to  A’n  Angstrom 
system  within  the  spectral  range  of  451  nm  to  725  nm.  The  laser  radiation  for  this  excitation 
is  derived  from  a  doubled  Nd:YAG  laser  which  pumps  a  pulsed  dye  laser.  The  output  of 
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the  dye  laser  is  doubled  and  mixed  with  the  residual  radiation  of  the  Nd:YAG  laser  at 
1064  nm. 

An  experiment  was  performed  on  a  CO-air  diffusion  flame  [9].  The  lowest  mole- 
fraction  measured  in  that  flame  was  14%  or  S^TO^m*3  at  1900  K.  Based  on  this  result, 
the  detection  limit  was  estimated  to  be  4.0,1020  m'3  at  300  K  [9],  This  exceeds  the  calcu¬ 
lated  concentration  of  CO  molecules  produced  in  a  single  pulse  [3].  Thus,  it  may  be  con¬ 
cluded  that  the  present  two-photon  LIF  technology  does  not  offer  the  sufficient  sensitivity 
required  for  the  detection  of  CO  traces  generated  from  atmospheric  COz  within  a  single¬ 
pulse  of  the  ArF  laser. 

As  an  alternative  to  the  conventional  Time-of-Flight  technique  we  proposed  [3]  to 
investigate  a  new  approach  in  which  the  tracer-generation  and  tracer-probing  pulses  con¬ 
verge  into  a  single  pulse.  In  this  approach  the  tracer  CO  molecules  are  generated  within 
the  high  intensity  ArF  laser  beam.  The  intense  radiation  may  further  dissociate  the 
molecule,  in  a  two-photon  process,  into  a  C(21D)  atom  and  an  0(3P)  atom.  The  ArF  laser 
at  A  =  193.1nm  is  resonant  with  the  transition  from  the  C(21D)  level  to  the  0(3^)  level. 
Therefore,  excitation  of  the  carbon  atoms  is  extremely  efficient  and  may  take  place,  due  to 
power-broadening,  at  wavelengths  which  are  slightly  off  resonance.  In  practice  most  of  the 
energy  of  a  broad-band  ArF  laser  is  coupled  into  this  transition  [5],  The  excited  atoms 
finally  radiate  at  247.8.  Although  this  laser-induced  fluorescence  is  a  four-photon  process  it 
benefits  from  the  high  transition-probability  in  the  carbon  atom.  Therefore,  at  tight  focus¬ 
ing  conditions,  the  intensity  of  the  radiation  at  A =247.8  nm  is  comparable  to  other  radia¬ 
tive  processes  induced  in  air  by  the  ArF  laser  [1]. 

The  intensity  of  the  fluorescent  radiation  in  this  four-photon  process  is  proportional 
to  the  density  of  CO  molecules.  Two  competing  mechanisms  determine  the  density  of  these 
molecules  along  the  laser  beam  path:  The  dominant  mechanism  is  the  generation  of  CO 
molecules  by  the  dissociation  of  C02.  However,  High  speed  air  flow  across  the  laser  beam 
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removes  some  of  the  generated  CO  molecules.  This  removal  process  which  competes  with 
the  generation  process  reduces  the  carbon-fluorescence  intensity.  Therefore,  the  reduced 
fluorescence  intensity  may  serve  for  the  measurement  of  the  air-flow  velocity  across  the 
beam. 

In  our  earlier  work  [3]  we  derived  the  rate  equations  which  describe  the  time- 
dependent  concentration  of  the  CO  molecules,  the  concentration  of  the  radiating  carbon 
atoms  and  the  concentration  of  the  carbon  atoms  which  were  trapped  in  a  metastable  state 
or  were  chemically  removed  from  the  flow.  Normally,  rate  equations  include  only  time- 
dependence.  They  do  not  include  any  spatial  dependence  nor  can  they  account  for  the  drift 
of  species  by  the  moving  air.  To  include  the  effect  of  gas  flow  across  the  laser  beam  the 
time-derivative  operator  d/dt  must  be  replaced  with  the  material-derivative  operator  D/Dt 
[10]: 

D/Dt  =  a/at  +  U-v  (2) 

where  U  is  the  flow-velocity  vector  and  V  is  the  gradient  operator.  The  velocimetry  tech¬ 
nique  described  here  is  sensitive  mostly  to  the  flow  rate  across  the  beam.  It  is  insensitive  to 
drift  along  the  laser  beam.  To  increase  sensitivity,  the  narrow  dimension  of  the  typically 
rectangular  ArF  beam,  will  be  aligned  with  the  cross-flow  (see  figure  2).  Therefore,  the 
material  derivative  operator  may  be  reduced  to  a  one  dimensional  operator: 

D/Dt=  a  fat  +  u  *a/a  (3) 

where  u  is  the  cross-flow  velocity. 

By  introducing  the  material-derivative  operator  into  the  rate  equations  they  become 
spatially  dependent.  Thus,  the  initial  conditions  and  all  the  variables,  such  as  the  laser- 
beam  intensity,  must  also  be  given  cheir  spatial  distribution.  Since  the  initial  density  of 
COz  is  homogeneous,  the  laser-beam  intensity  is  the  o:dy  forcing  term  with  a  non-uniform 
spatial  distribution.  A  typical  ArF  laser  beam  is  emitted  with  a  cross-section  of  1x2.5cm  and 
has  an  approximately  uniform  intensity  distribution  across  both  dimensions.  However, 
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when  a  plane-wave  radiation  with  an  energy  Eo  and  with  a  uniform  intensity  distribution  is 
focused,  the  intensity  distribution  at  the  focal  plane  is  [11]: 
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where  x  an  y  are  the  coordinates  in  the  focal  plane,  AX  and  AY  are  the  unfocussed 
beam  width  and  height  respectively,  f  is  the  focusing  lens  focal  length,  r  is  the  pulse  dura¬ 


tion  and: 


.  .  .  sinjrx 

smc(x)  = - 

tX 


In  most  applications  to  be  described,  the  x  dimension  will  be  much  narrower  than  the  y 
dimension.  Hence,  the  intensity  distribution  may  be  simplified  by  describing  the  intensity 
distribution  in  the  y  direction  as  being  uniform  and  assigning  for  the  intensity  distribution 
in  the  x  direction  the  following  equation: 


T  ■  Eq*X 

r*AYAf 


sine2 


To  further  simplify  the  analysis,  the  distribution  in  the  x  direction  may  be  replaced  by  a 
Gaussian  intensity-distribution.  After  normalization,  the  intensity  distribution  is: 

exp(iS>  (7) 

where  AX  is  the  coordinate  in  the  x  direction  at  which  the  intensity  falls  to  1/e2  of  its  cen¬ 
terline  value  and  I0  is  the  intensity  at  the  beam  centerline.  Usually  2aX  is  considered  as 
the  beam  width.  Introducing  also  the  time  dependence,  which  was  assumed  to  be  repre¬ 
sented  by  an  exponential  decay  [3],  we  obtain: 

1  “  exp(f->  (8) 


Here  E0  is  the  pulse  energy  and  r  is  the  time  at  which  the  intensity  decays  to  1/e  of  its  ini¬ 
tial  value.  Usually,  t  is  considered  as  the  pulse  duration.  The  rate  equations  [3],  together 
v/ith  the  material-derivative  operator  (equation  (3)  and  equation  (8))  describe  the  time- 
dependent  spatial-distribution  of  the  dissociation  products  and  of  the  flux  of  fluorescence 
photons  at  247.8  nm  produced  by  a  focused  ArF  laser  beam  which  is  probing  the  air-flow. 

2.2  RESULTS 

2.2.1  FLUCTUATIONS  IN  THE  LASER  INTENSITY  PROFILE 

Equation  (8)  presents  a  realistic  distribution  of  the  energy  across  the  laser  beam 
profile.  This  equation  is  sufficient  for  accurately  predicting  the  total  radiation  emitted  in 
linear  fluorescence  processes.  However,  in  interactions  which  depend  non-linearly  on  the 
laser  intensity,  even  slight  variations  in  the  intensity  distribution,  may  create  large  fluctua¬ 
tions  in  the  induced  fluorescence  energy.  Like  most  gas  lasers  the  intensity  profile  of  the 
ArF  laser  may  present  some  pulse  to  pulse  fluctuations  which  are  induced  by  irregularities 
in  the  excitation  discharge.  To  estimate  the  extent  of  this  effect  we  introduced  two  alterna¬ 
tive  intensity  distribution  functions  and  calculated  the  total  carbon  fluorescence  which  will 
be  induced  by  laser  beams  with  such  distributions  while  all  other  parameters  remain 
unchanged. 

To  simulate  the  effect  of  ’hot  spots’  we  presented  the  laser  beam  as  a  superposition 
of  two  Gaussian  distributions  offset  from  the  beam  centerline.  The  new  distribution  func¬ 
tion  was  normalized  such  that  the  total  energy,  E0,  was  equal  to  the  total  energy  delivered 
by  the  single  Gaussian  distribution  beam.  The  new  distribution  which  presents  two  intensity 
spikes  is: 
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(9) 


where  Xq  is  the  displacement  of  the  Gaussian  distribution  from  the  beam  centerline.  In  ad¬ 
dition,  we  also  simulated  a  homogeneous  energy  distribution  in  which  all  the  energy  is 
uniformly  distributed  within  a  spot  of  dimensions  of  2AX2AY.  The  distribution  function 
for  this  configuration  is: 


i  =  E° 


2nX.AY7eXP(f) 


(10) 


Also  here  the  total  energy  under  the  distribution  curve  is  E  . 

A  numerical  code  was  developed  [3]  and  used  to  solve  the  time  dependent  and  spa¬ 
tially  varying  rate  equations.  In  this  code,  each  solution  step  included:  (a)  an  explicit 
forward-difference  solution  [12]  of  one  half-step  in  the  convection  term  while  the  time  is 
held  constant,  followed  by  (b)  a  predictor-corrector  solution  of  one  time-step  while  the  the 
convection  term  is  frozen,  followed  by  (c)  the  remaining  half  step  in  the  explicit-forward- 
difference  convection  term. 

The  stability  of  this  algorithm  is  limited  by  the  explicit-forward-difference  steps,  with 
a  stability  criterion  of  [12]: 

u'6t/£x  <  <  1  (10) 

where  sx  and  St  are  the  sizes  of  the  spatial  and  time  steps  in  the  algorithm. 

The  integration  limits  in  the  calculation  were  |x|  <4aX  and  t<4r.  Beyond  these 
limits,  the  intensity  of  the  laser,  for  all  distribution  functions,  was  less  then  1.8%  of  its  maxi¬ 
mum.  The  number  of  steps  in  the  x  direction  was  set  to  200  while  the  number  of  time  step 
was  selected  individually  for  each  velocity  using  equation  (10).  Owing  to  the  uniform  inten- 
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sity  distribution  in  the  y  direction,  the  beam  height,  AY,  could  be  selected  arbitrarily.  It 
was  chosen,  however,  such  that  E0= lOOmJ  for  a  prescribed  intensity. 

To  illustrate  the  effect  of  an  air-flow  across  the  laser  beam  on  the  density  of  the 
various  carbon  species  and  on  the  total  fluorescence  intensity  at  247.8  nm  we  solved  the  al¬ 
gorithm  for  a  case  in  which  Efl=100  mJ,  r=20nsec,  AX  =10^m,  AY  =  .399  mm,  and  for 
velocities  ranging  from  u  =  0  to  2500  m/sec.  This  focusing  configuration  corresponds  ap¬ 
proximately  to  the  highest  photon-flux  which  may  be  obtained  without  inducing  air  break¬ 
down  [1,2].  The  peak  intensity  for  the  single  Gaussian  distribution  was  I  =  T1015  watts/m2. 
It  was  assumed  in  the  calculation  that  the  narrow  dimension,  i,e.  AX,  is  parallel  to  the 
main  flow-velocity  component. 

The  photon-flux  emitted  by  the  excited  C  atoms  when  v=0  m/sec  and  when  the  in¬ 
tensity  distribution  is  described  by  a  single  Gaussian  is  1.435-1013  photons/m2.  We  assume 
that  the  radiation  is  collected  with  a  1mm  resolution.  We  further  assume  that  the  collection 
of  the  light  is  accomplished  by  an  f/2  lens,  that  the  fluorescence  is  separated  from  the  back¬ 
ground  radiation  by  a  band-pass  filter  with  an  efficiency  of  r?=0.5  and  that  the  photo¬ 
detector  quantum  efficiency  is  »?d=0.15.  With  these,  the  number  of  photo-electrons  ob¬ 
served  at  A=247.8nm  by  focusing  lOOmJ  of  radiation  at  193. lnm  into  atmospheric  air  is 
N  =16871  photons.  Similarly  when  u=200  m/sec  Npe=  15490  photons  and  when  u=2400 
m/sec  =4546  photons. 

Figure  3  describes  the  number  of  photo-electrons  detected  at  247.8  nm  vs  the  flow 
velocity  for  the  three  ditfereut  intensity  distribution  functions.  For  all  distributions 
AX=10/xm.  The  solid  line  describes  the  total  fluorescence  observed  when  the  intensity 
distribution  is  a  single  Gaussian,  the  doted  line  corresponds  to  the  double  Gaussian  dis¬ 
tribution  where  x0=5/un  while  the  dashed  line  corresponds  to  the  uniform  distribution.  To 
simplify  the  comparison,  hi  curves  were  normalized  by  16871  photo-electrons  which  is  the 
count  when  u=0  and  the  energy  distribution  is  described  by  a  single  Gaussian. 
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RELATIVE  INTENSITY 


T  -  20  nsec.,E0  =  100  mJ,  ( u=0 )  =  16871  count 


VELOCITY  m  /sec 

Figure  3.  Total  fluorescence  at  247.8  nm  vs  the  flow  velocity  for  the  three  different  in¬ 
tensity  distribution  functions.  The  solid  line  describes  the  total  fluorescence 
observed  when  the  intensity  distribution  is  a  single  Gaussian,  the  dotted  line 
corresponds  to  the  double  Gaussian  distribution  where  xf)=5nm  and  the 
dashed  line  corresponds  to  the  uniform  distribution.  For  all  distributions 
AX  =  10/nn. 
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It  is  seen  from  this  figure  that  when  the  distribution  of  the  intensity  across  the  laser 
beam  varies,  the  number  of  observed  photons  will  vary  as  well.  It  is  impossible  to  predict 
the  intensity  distribution  patterns  when  the  discharge  is  unstable.  Therefore,  we  assumed 
that  any  one  of  the  distributions  which  are  described  here  may  occur  at  random.  In  such  a 
case  the  fluorescence  intensity  will  fluctuate  randomly  as  well.  At  low  velocities  the 
measurement  error  may  be  as  high  as  100%.  When  the  air  velocity  is  1000  m/sec  the  num¬ 
ber  of  detected  photons  induced  by  the  beam  with  the  single  Gaussian  distribution  equals 
to  the  number  of  photons  which  are  induced  by  the  beam  with  the  uniform  distribution 
when  the  velocity  is  604  m/sec  and  to  the  number  of  photons  which  are  induced  by  the 
beam  with  the  double  Gaussian  distribution  when  the  velocity  is  875  m/sec.  This  cor¬ 
responds  to  an  uncertainty  in  the  velocity  measurement  which  may  be  as  high  as  40%. 
Similar  absolute  errors  will  be  observed  when  the  air  speed  is  2500  m/sec.  Here  however 
the  relative  uncertainty  is  less  than  16%. 

2.2.2  THE  EFFECT  OF  TEMPERATURE  AND  DENSITY  ON  THE  MEASURE¬ 
MENT  ACCURACY 

Most  laser  induced  fluorescence  interactions  are  temperature  dependent.  In  equi¬ 
librium,  the  population  of  individual  states  vary  with  temperature  according  to  the 
Boltzmann’s  distribution.  Therefore,  when  individual  states  are  excited  by  a  narrowband 
laser  the  total  fluorescence  reflects  the  temperature  dependence  of  the  population  of  the 
initial  state.  In  addition,  at  elevated  temperatures,  dissociation  of  the  ground  state  and  of 
the  excited  molecules  reduces  the  density  of  the  interacting  molecules  and  thereby  the  in¬ 
tensity  of  the  observed  fluorescence.  Finally,  Doppler  and  pressure  broadening  are  tem¬ 
perature  dependent  and  may  affect  the  coupling  between  a  narrowband  laser  and  the 


14 


probed  molecules.  We  investigated  the  effect  of  temperature  on  the  fluorescence  intensity 
emitted  by  C  atoms  following  a  four-photon  dissociation  and  excitation  of  C02. 

The  dissociation  of  COz  and  CO  by  the  ArF  laser  is  broadband  in  nature.  The  in¬ 
teraction  does  not  depend  on  the  excitation  of  an  individual  line.  Instead  it  is  represented 
by  a  transition  from  the  ground  state  into  an  upper  sb  te  continuum.  Therefore,  the  dis¬ 
sociation  process  which  does  not  depend  on  the  population  of  an  individual  level  in  the 
ground  state  is  independent  of  temperature.  The  interaction  with  the  atomic  carbon  is  line 
specific.  Although  this  is  expected  to  be  an  atomic  narrowband  transition,  it  appears,  due  to 
power  broadening,  almost  as  broad  as  the  tuning  range  of  the  ArF  laser  [3].  Therefore,  the 
effects  of  the  Doppler  and  pressure  broadening  mechanisms  is  expected  to  be  negligible.  In 
addition,  all  carbon  atoms  here  are  prepared  at  the  2^  state  by  the  dissociation  process. 
Therefore,  the  population  density  at  this  state  is  determined  by  the  laser  parameters  and  by 
the  concentration  of  the  COz  molecules  and  not  by  the  equilibrium  temperature.  Further¬ 
more,  most  of  these  atoms  are  consumed  by  various  recombination  processes  which  occur 
within  one  collision.  As  a  result,  the  population  of  this  state  cannot  be  equilibrated  and  it 
remains  temperature  independent. 

It  is  expected  that  only  dissociation  of  C02  or  CO  which  occurs  at  elevated  tempera¬ 
tures  can  affect  the  total  fluorescence  emitted  by  the  C  atoms.  This  effect  is  likely  to  be 
small  at  room  conditions  but  may  become  significant  in  high  temperature  stagnation  points 
or  in  reacting  flows. 

To  examine  this  effect  we  calculated  the  chemical  equilibrium  of  dry  air  using  the 
STANJAN  iterative  computer  program  which  was  developed  by  W.C.Reynolds  from  Stan¬ 
ford  University. 

The  program  models  the  gas  as  a  mixture  of  ideal  gases  and  possibly  condensed 
phases  in  the  form  of  liquid  droplets  or  solid  particulates.  The  thermodynamical  data  base 
is  taken  from  the  JANAF  tables  and  is  structured  accordingly.  The  user  is  expected  to  iden- 
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tify  the  species  of  the  mixture,  the  relative  concentration  of  the  atomic  components,  and 
two  independent  thermodynamic  properties  of  the  state  at  which  the  equilibrium  is  sought 
(e.g.  temperature  and  pressure).  The  calculation  of  the  equilibrium  is  based  on  the  method 
of  element  potentials  in  which  exact  equations  for  the  gas-phase  mole  fractions  are  derived 
in  terms  of  Lagrange  multipliers  and  the  total  number  of  moles  is  adjusted  to  meet  the 
various  atomic  concentrations  and  to  render  the  sum  of  the  mole  fractions  to  unity.  Thus,  in 
mixtures  which  involve  N  species  but  only  M  atoms  (where  N  >  >  M),  the  program  requires 
the  adjustment  of  only  M+l  parameters.  This  is  in  contrast  to  most  equilibrium  analysis 
techniques  where  all  the  unknown  mole  fractions  need  to  be  adjusted. 

We  used  a  modified  IBM  PC  version  of  this  program.  We  applied  it  for  the  equi¬ 
librium  analysis  of  air  which  included  the  following  constituents  with  their  initial  molar 
fraction:  0.7809  N2,  0.2095  Oz,  0,0093  Ar,  and  0.0003C02.  The  program  was  tested  by  com¬ 
paring  the  calculated  equilibrium  mixture  of  air  at  1  atm  and  at  temperatures  of  1500  K 
and  2500  K  with  the  equilibrium  mixture  obtained  by  using  the  CEC72  equilibrium 
program  which  was  developed  by  Svehla  and  McBride  from  NASA  Lewis.  The  results  of  of 
the  CEC72  program  agreed  with  the  results  of  the  STANJAN  program  within  1%.  This 
confirmed  the  accuracy  of  the  STANJAN  program. 

The  variation  of  the  molar  fraction  of  CO  with  temperature  is  presented  in  figure  4 
while  the  temperature  dependence  of  C02  is  presented  in  figure  5.  The  pressure  for  these 
calculations  was  kept  at  1  atmosphere  while  the  temperature  was  varied  between  1000  - 
3000  K,  The  relative  concentration  of  CO  is  below  10‘6  throughout  the  whole  temperature 
range  of  interest.  This  concentration  is  small  relative  to  the  concentration  of  CO  following 
laser  induced  C02  dissociation.  The  molar  fraction  of  C02  is  shown  in  figure  5  to  be  inde¬ 
pendent  of  temperature  when  the  air  temperature  is  below  2000  K.  At  higher  tempera¬ 
tures,  dissociation  takes  place,  the  concentration  of  COz  decreases  while  the  concentration 
of  C  increases.  Therefore,  below  2000  K,  when  the  COz  concentration  is  temperature 
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The  variation  with  temperature  of  the  molar  fraction  of  CO  in  air  in  terms  of 
<0gw  (10’3-[COJ). 


17 


o  |  t  i  i  r~ 

1000  1400  1800  2200  2600  3000 

TEMPERATURE -K 


Figure  5:  The  variation  with  temperature  of  the  molar  fraction  of  C02  in  air  in  terms 

of  lO^COJ. 
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independent,  the  fluorescence  intensity  emitted  by  C  atoms  at  247.8  nm  is  expected  to  be 
temperature  independent.  Therefore,  a  single  measurement  of  the  total  fluorescence  at 
247.8  nm  and  of  the  gas  density  will  yield  the  flow  velocity  across  the  laser  beam.  Note  that 
both  the  C02  and  CO  concentrations  vary  linearly  with  the  initial  air  density.  Therefore, 
these  results  are  expected  to  be  valid  also  at  other  densities.  We  repeated  our  calculations 
for  P=0.01  atm  and  for  P  =  10  atm  and  found  that  also  at  these  pressures  the  COz  con¬ 
centration  was  independent  of  temperatures  when  the  air  temperature  is  below  2000K. 

2.2.3  EXPERIMENTAL  TESTING  OF  THE  FEASIBILITY  OF  THE  TECHNIQUE 

The  single  pulse  carbon  velocimetry  is  expected  to  be  sensitive  to  the  velocity  of 
high  speed  flows.  We  tested  the  feasibility  of  this  technique  by  comparing  the  total  fluores¬ 
cence  emitted  by  the  dissociation  products  of  atmospheric  C02  in  stationary  air  with  the  to¬ 
tal  fluorescence  emitted  in  supersonic  jets. 

The  experimental  setup  included  a  broadband  ArF  laser  which  emits  250  mJ/pulse 
at  A  =  193  nm  at  a  rate  of  10Hz.  The  beam  was  focused  by  a  lens  with  a  focal  length  of  150 
mm.  The  spot  size  was  estimated,  using  a  Polaroid  film  burning,  to  be  less  than  100  nm. 
This  was  the  tightest  focusing,  using  this  amount  of  energy,  which  could  be  obtained 
without  inducing  an  air  breakdown.  The  induced  fluorescence  was  collected  by  an  f/2  lens 
and  was  imaged  on  the  slits  of  a  Jarrel-Ash  0.5  m  monochromator  which  was  equipped  with 
a  1180  grooves/mm  grating  blazed  at  5000  A  and  a  solar  blind  photomultiplier  tube 
(PMT).  The  output  of  the  PMT  was  monitored  by  a  100  MHz  oscilloscope.  The 
monochromator  slits  were  set  at  250  Mm.  This  was  the  largest  slit  width  which  still  allowed 
to  resolve  the  247.8  nm  carbon  line  from  the  adjacent  oxygen  fluorescence  line. 

To  create  a  supersonic  jet  we  designed  a  1/8"  orifice  which  could  be  mounted  on  the 
end  of  a  tube  which  was  connected  to  a  high  pressure  air  cylinder.  The  expansion  ratio  was 
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designed  such  that  for  stagnation  conditions  of  150  psi  and  25  C  and  for  an  isentropic  fully 
expended  flow  the  exit  Mach  number  was  M=2.2.  To  assure  that  the  jet  is  correctly  aligned 
both  with  the  tightly  focused  beam  and  with  the  collection  cone  we  used  initially  a  helium 
flow.  When  the  jet  was  aligned  correctly  the  signal  both  due  to  the  carbon  line  and  due  to 
the  oxygen  fluorescence  diminished.  When  the  supply  line  was  reconnected  to  the  air 
supply  the  signal  recovered. 

We  performed  two  measurements  of  the  fluorescence  intensity:  a  measurement  with 
the  laser  beam  propagating  through  the  stationary  room  air  and  a  measurement  with  the 
beam  propagating  through  the  supersonic  jet.  The  intensity  measurements  in  both  experi¬ 
ments  were  approximately  identical  although  the  high  speed  flow  across  the  beam  wt  3  ex¬ 
pected  to  reduce  the  fluorescence  intensity.  Changing  the  stagnation  pressure  and  hence 
the  total  density  also  did  not  have  any  noticeable  effect  on  the  observed  fluorescence. 

In  a  subsequent  work  under  our  NASA  contract  we  studied  some  of  the  properties 
of  the  non-linear  process  in  C02.  It  was  concluded  that  when  the  laser  fluence  exceeds  1.5 
J/cm2,  ionization  and  quenching  effects  interfere  with  the  dissociation  and  excitation 
processes.  As  a  result,  the  otherwise  third  order  process  may  appear  as  almost  independent 
of  laser  intensity  or  gas  density.  The  maximum  fluence  in  the  present  experiments  was  2000 
J/cm2.  This  fluence  exceeds  the  minimum  fluence  of  1.5/ J/cm2  at  which  the  non-linear 
process  follows  the  predictable  third-order  behavior.  On  the  other  hand,  the  minimum 
fluence  which  is  required  to  obtain  an  accuracy  of  at  least  5%  is  66  J/cm2  [3].  The  dif¬ 
ference  between  the  minimum  fluence  required  for  high  accuracy'  and  the  maximum 
fluence  at  which  ionization  and  quenching  effects  are  negligible  is  too  large  and  can  not  be 
reconciled. 

It  is  advantageous  to  use  high  fluence  for  these  measurements.  However,  the  inter¬ 
pretation  of  the  fluorescence  measurements  for  single  pulse  carbon  velocimetry  must  in¬ 
clude  also  an  analysis  of  the  interfering  ionization  and  quenching  effects.  This  analysis  is 
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more  complex  than  the  interpretation  of  the  results  of  the  third  order  process.  In  addition, 
the  normalization  of  signals  can  not  rely  on  a  simple  calibration  cell.  Therefore,  we  believe 
that  the  ionization  and  quenching  effects  associated  with  the  required  high  fluence  make 
the  single  pulse  carbon  velocimetry  too  complex  for  most  wind  tunnel  applications. 
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An  alternative  to  the  CO  velocimetry  is  a  velocity  measurement  technique  which 
relies  on  the  Doppler  shift  experienced  by  molecules  moving  in  the  gas  stream.  This  tech¬ 
nique  was  originally  proposed  and  demonstrated  by  one  of  the  co-workers  of  this  project 
(Zimmermann  [13]).  In  this  technique,  the  molecules  are  excited  by  a  tunable  dye-laser  and 
their  excitation  is  marked  by  the  fluorescence  decay  from  the  excited  state.  The  difference 
between  the  frequency  needed  to  excite  stationary  molecules  and  the  frequency  at  which 
the  moving  molecules  are  excited  is  the  Doppler  shift.  This  shift  is  proportional  to  the  gas 
velocity. 

Several  schemes  for  the  measurement  of  the  Doppler  shift  have  been  proposed. 
Zimmermann  used  successive  counter-propagating  beams  to  measure  the  frequency  dif¬ 
ference  between  the  up-shifted  and  down-shifted  frequencies  of  the  D  lines  of  Na  atoms 
seeded  into  a  He  flow  [13]. 

She  [14],  Exton  [15]  and  others  proposed  to  replace  the  LIF  technique,  which  was 
used  to  measure  the  Doppler-shift,  with  a  coherent  Raman  process  such  as  the  InvCiSe- 
Raman-Spectroscopy  (IRS)  or  CARS.  Both  these  coherent  Raman  schemes  require  a  fixed 
frequency  laser  and  a  synchronously  pumped  narrow-band  tunable  laser  which  is  tuned  to 
excite  the  Raman  process.  When  the  two  laser  beams  which  are  required  for  the  excitation 
of  the  Raman  process  are  counter  propagating,  the  Doppler  shift  sensed  by  the  molecules 
and  the  measurement  sensitivity  are  doubled.  Co-propagating  beams  are  usually  used  in 
stimulated  Raman  spectroscopy  to  suppress  the  Doppler  shifts.  Both  the  LIF  and  the 
Raman  Doppler  techniques  are  slow.  The  time  required  to  tune  the  probe  laser  until  the 
Doopler  shift  is  measured  may  last  up  to  few  minutes  per  point  and  therefore  is  suitable 
only  for  time-averaged  measurements  [13], 
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Hiller  et  al.[16]  proposed  a  method  which  removes  this  difficulty.  In  their  approach, 
two  counter-propagating  laser  beams  are  tuned  to  the  high  frequency  wing  of  the  absorp¬ 
tion  line  of  I2  molecules  seeded  in  a  N2  subsonic  flow.  In  the  absence  of  flow,  both  beams 
enco.ter  the  same  absorption  cross-section  of  the  I2  molecules  and  induce  the  same 
fluorescence  intensity.  However,  when  the  absorbing  molecules  move  along  the  laser 
beams  together  with  the  gas  flow,  the  absorbed  radiation  is  Doppler  shifted:  the  frequency 
of  the  radiation  of  the  beam  which  propagates  opposite  to  the  flow  is  Doppler  up-shifted 
and  is  therefore  further  away  from  the  absorption  line-center.  Therefore,  fluorescence  in¬ 
tensity  induced  in  molecules  which  move  opposite  to  the  beam  decreases.  The  other  laser 
beam  -  which  propagates  in  the  direction  of  the  N2  flow  -  is  experienced  by  the  moving 
molecules  as  Doppler  down-shifted.  Therefore,  the  fluorescence  induced  by  this  beam  will 
be  more  intense.  The  difference  in  intensities  between  these  two  LIF  processes  has  been 
shown  by  Hassa  et  al.  [17]  to  relate  to  the  flow  velocity. 

To  further  simplify  this  Fluorescence  Difference  (FD)  velocimetry  technique  and  to 
eliminate  the  need  for  seeding  we  proposed  [3]  to  use  an  absorption  line  in  the  Schumann- 
Runge  band  of  Oz  for  monitoring  Doppler  shifts  experienced  by  these  molecules  in  air 
flows.  This  technique  benefits  from  the  coincidence  between  the  frequency  of  the  ArF  laser 
and  the  absorption  frequency  of  Oz  in  the  Schumann-Runge  band.  In  this  transition  the  Oz 
molecule  is  excited  from  the  X3EQ'  electronic  ground  state  to  the  B3SU*  electronic  ex¬ 
cited  state  [18].  This  approach  is  compatible  with  our  ongoing  effort  of  developing  laser 
techniques  for  the  measurement  of  air  temperature  and  density  in  hypersonic  wind  tunnel 
flows. 

The  measurement  uncertainty  of  all  Doppler  techniques  is  determined  primarily  by 
the  spectral  line  width  of  the  laser  and  by  the  width  of  the  transition  line  of  the  probed 
molecule.  The  reported  line-width  of  the  only  commercially  available  tunable  ArF  laser  is 
A.L  =  9  GHz  [19],  The  natural  line-width  Ai/u  of  the  v"(0)  to  v’(4)  transition  of  02, 
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which  may  be  excited  by  the  ArF  laser,  is  125.4  GHz  [20].  This  extremely  large  width  is  due 
to  the  predissociation  of  the  upper  state.  For  comparison,  the  Doppler  shift  A^D  of  the 
moving  flow  is  [21]: 


Av0  =  vu/c  =u/A  (11) 

For  a  flow  of  u  =  1000  m/sec  and  for  a  wavelength  of  A  =  193.3  nm,  the  Doppler  shift  is 
5.17  GHz.  This  shift  is  small  relative  to  the  v"(0)  to  v*(4)  transition  line-width  and  is  com¬ 
parable  to  the  laser  line-width.  Although  the  application  of  an  ArF  laser  beam  for  the 
measurement  of  the  Doppler  shift  experienced  by  the  v"(0)  to  v’(4)  transition  of  02  may 
yield  a  qualitative  velocity  measurement  it  is  unlikely  to  result  in  an  accurate  quantitative 
measurement.  The  accuracy  of  the  measurement  may  be  enhanced  by  selecting  a  laser 
wavelength  which  is  resonant  with  an  Oz  transition  which  presents  a  narrower  line-width. 
The  predissociation  life-time  of  the  v’(13)  K=13  level  in  the  B3£L  electronic  state  cor¬ 
responds  to  a  natural  line-width  of  4.58GHz  [20].  Transitions  from  the  ground  state  to  this 
state  result  in  a  line  which  is  among  the  narrowest  in  this  manifold  of  vibrational  states  and 
therefore  is  the  most  suitable  for  the  Doppler  velocimetry  technique.  This  line  is  resonant 
with  the  wavelength  of  the  ArF  laser  which  was  Raman  shifted  with  a  H2  Raman  cell  [22]. 
This  device,  which  is  available  commercially,  uses  high  pressure  H2  as  the  Raman  gain 
medium  to  shift  the  wavelength  of  the  beam  by  4155  cm'1.  If  the  ArF  laser  beam  is  tuned  to 
193.138  nm  the  wavelength  of  the  Raman  shifted  beam  will  be  A  =  178.790  nm  which  is 
resonant  with  the  v"(0)  to  v’(13)  R13  transition. 

In  a  previous  work  [3],  we  calculated  analytically  the  measurement  uncertainty 
which  is  associated  with  this  velocimetry  technique.  The  calculation  assumed  interaction 
w'th  only  one  rotational  transition  in  the  v’(O)  to  ^(13)  band.  However,  the  Raman  shifted 
laser  beam  is  also  resonant  with  the  vh(0)  to  v*=  14,15,16  and  17  transitions.  Some  of  these 
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transitions  overlap  the  v"(0)  to  v’(13)  band.  Therefore,  an  accurate  prediction  of  the  tem¬ 
perature  dependence  of  this  measurement  requires  the  inclusion  of  all  the  transitions  which 
are  resonant  with  the  laser.  We  have  modified  a  model  which  was  developed  by  us  under  a 
separate  contract  [2].  The  modified  model  describes  the  LiF  process  in  molecular  oxygen 
excited  by  a  Raman  shifted  beam.  The  model  calculates  the  intensity  of  the  total  fluores¬ 
cence  emitted  by  an  air  stream  moving  along  and  opposite  to  a  Raman  shifted  ArF  laser 
beam.  Both  the  air  velocity,  temperature  and  density  can  be  included  as  parameters.  In  the 
following  sections  we  describe  the  model  and  the  results  of  our  calculations. 

3.1  THE  MODEL  OF  LIF  IN  OXYGEN 

A  computer  model  which  describes  the  absorption  of  an  ArF  laser  beam  by  air  and 
the  subsequent  induced  fluorescence  has  been  developed.  The  code  predicts  the  effects  of 
absorption  by  room  air  of  a  Raman  shifted  ArF  laser  beam  with  a  prescribed  initial 
spectral  distribution  as  it  propagates  from  the  laser  to  the  sample  volume.  The  attenuation 
at  each  wavelength  in  the  laser  bandwidth  is  evaluated  by  adding  the  contributions  of  all 
participating  transitions  that  originate  from  the  ground  electronic  state.  The  total  fluores¬ 
cence  energy  which  is  emitted  by  the  gas  in  the  observed  sample  volume  is  then  calculated 
using  the  spectral  distribution  of  the  incident  beam  at  the  sampling  point. 

3.1.1  The  Spectroscopic  Features  of  LIF  in  Oxygen 

The  calculation  of  absorption  includes  all  significant  rotational  transitioas 
which  fall  within  the  spectral  tuning  range  of  the  ArF  laser,  including  the  following  vibra¬ 
tional  bands:  from  v"(0)  to  v’(13)  v’(14),  v’(15),  ^(16),  and  v’(17).  No  other  transitions  from 
v"(0),  v"(l)  and  v"(2)  were  found  to  be  resonant  with  the  Raman  shifted  ArF  laser 
throughout  its  tuning  range.  The  term  values  of  the  ro-vibronic  levels  in  the  B3£'u  upper 
electronic  state  were  calculated  using  the  equations  and  parameters  of  Cheung  et  al  [23] 
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which  include  the  fine-splitting  of  the  rotational  states.  Those  upper-state  term  values, 
together  with  the  term  values  of  the  vH(0)  and  vM(l)  levels  of  the  X3E*9  lower  electronic 
state  which  are  given  by  Veseth  and  Lofthus  [24]  reproduce  the  measured  transition  fre¬ 
quencies  [25]  for  all  tra  ditions  originating  from  the  v"(0)  and  v"(i)  state.  For  transitions 
from  the  v"(2)  which  were  not  given  by  Ref.  24,  the  expressions  of  Endo  and  Mizushima 
[26]  for  the  dependence  of  the  parameters  on  vibrational  quantum  number  were  used. 
However,  their  constants  were  adjusted  to  reproduce  the  parameters  of  the  v"(0)  and  v"(l) 
states  reported  by  Ref.  24.  Thus,  for  the  rotational  term  values  of  all  the  vibrational  states 
in  the  X3E"g  electronic  state,  the  following  expressions  were  used: 

Bv=  1.445584  -0.015778(v+ 0.5) 

Dv  =  4.8425xl0'6  +  3.3xl0”a(v  +  0.5) 

A  =  1.9855  +  7x10'4(v+0.5)  (12) 

7  =  -8.525xl0"3  -  7x10'5(v+0.5) 
ad  =  1.839x10-6 
70  =  -4.01xl0's 

where  Bv,  Dy,A,  7,  AD  and  70  are  all  specified  in  cm'1  and  are  the  same  notation  which  is 
used  in  Refs.  23-26.  The  calculated  wavelengths  for  the  transitions  from  v"(0)  and  v"(l) 
matched  the  measured  wavelengths  [25]  within  0.5  cm'1.  No  measurements  were  found  to 
compare  with  computed  results  for  transitions  originating  from  v"(2)  and  they  are  not  ex¬ 
pected  to  have  the  same  accuracy  as  the  results  for  transitions  from  the  lower  states. 
However,  since  only  the  transitions  from  v’(0)  were  found  to  be  resonant  with  the  Raman 
shifted  beam  this  will  not  affect  the  accuracy  of  this  calculation.  Table  1  lists  all  the  transi¬ 
tions  and  their  wavelengths  which  are  resonant  with  the  Raman  shifted  ArF  laser  beam. 
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FUNDAMENTAL 

RAMAN 

WAVELENGTH 

SHIFTED 

I 

A! 

W,L(A‘) 

FREQUENCY 

1 

1929.739 

1786.497 

55975.469 

2 

1929.681 

1786.447 

55977.039 

3 

1929.657 

1786.426 

55977.680 

1 

1931.378 

1787.901 

55931.516 

2 

1931314 

1787.846 

55933215 

3 

1931.286 

1787.822 

55933.980 

1 

1933.297 

1789.545 

55880.121 

2 

1933.228 

1789.486 

55881.965 

3 

1933.194 

1789.457 

55882.863 

1 

1935.499 

1791.432 

55821.266 

2 

1935.425 

1791.369 

55823.242 

3 

1935.386 

1791.335 

55824.277 

1 

1929.167 

1786.006 

55990.852 

2 

1929.051 

1785.907 

55993.949 

1 

1932.069 

1788.493 

55912.988 

2 

1931.946 

1788.388 

55916.285 

3 

1931.879 

1788.330 

55918.078 

1 

1935.281 

1791.245 

55827.082 

2 

1935.150 

1791.133 

55830.590 

3 

1935.075 

1791.069 

55832.590 

1 

1932.246 

1788.645 

55908.234 

2 

1932.041 

1788.469 

55913.742 

3 

1931.926 

1788.371 

55916.809 

1 

1936.284 

1792.104 

55800.328 

2 

1936.066 

1791.917 

55806.141 

3 

1935.939 

1791.809 

55809.508 

1 

1931.602 

1788.093 

55925.508 

2 

1931.270 

1787.809 

55934.391 

3 

1931.088 

1787.652 

55939.289 

1 

1936.215 

1792.046 

55802.152 

2 

1935.865 

1791.746 

55811.488 

3 

1935.665 

1791.574 

55816.836 

1 

1933.158 

1789.426 

55883.828 

2 

1932.644 

1788.986 

55897.598 

3 

1932.354 

1788.737 

55905.359 
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THE  P  BRANCH 


FUNDAMENTAL  RAMAN 


Y1 

JC 

K 

I 

WAVELENGTH 
A : 

SHIFTED 

W.L.(A*1 

FREQUENCY 

0 

13 

11 

1 

1930.625 

1787.256 

55951.703 

0 

13 

11 

2 

1930.572 

1787.210 

55953.125 

0 

13 

11 

3 

1930.553 

1787.194 

55953.625 

0 

13 

13 

■4 

1932.416 

1788.791 

55903.684 

0 

13 

13 

2 

1932.359 

1788.741 

55905.230 

0 

13 

13 

3 

1932.335 

1788.722 

55905.848 

0 

13 

15 

1 

A 

1934.488 

1790.566 

55848.266 

0 

13 

15 

2 

1934.425 

1790.512 

55849.949 

0 

13 

15 

3 

1934.397 

1790.488 

55850.695 

0 

13 

17 

1 

1936.842 

1792.583 

55785.430 

0 

13 

17 

2 

1936.774 

1792.524 

55787.250 

0 

13 

17 

3 

1936.741 

1792.496 

55788.129 

0 

14 

19 

1 

1930,516 

1787.162 

55954.621 

0 

14 

19 

2 

1930.409 

1787.071 

55957.500 

0 

14 

19 

3 

1930.358 

1787.027 

55958.867 

0 

14 

21 

1 

1933.551 

1789.763 

55873.316 

0 

14 

21 

2 

1933.436 

1789.665 

55876.391 

0 

14 

21 

3 

1933.377 

1789.614 

55877.961 

0 

14 

23 

1 

1936.894 

1792.627 

55784.055 

0 

14 

23 

2 

1936.771 

1792.521 

55787.336 

0 

14 

23 

3 

1936.704 

1792.464 

55789.109 

0 

15 

23 

1 

1929.976 

1786.699 

55969.121 

0 

15 

23 

2 

1929.793 

1786.543 

55974.016 

0 

15 

23 

3 

1929.702 

1786.465 

55976.469 

0 

15 

25 

1 

1933.782 

1789.961 

55867.137 

0 

15 

25 

2 

1933.588 

1789.795 

55872.328 

0 

15 

25 

3 

1933.485 

1789.707 

55875.078 

0 

16 

27 

1 

1933.008 

1789.298 

55887.836 

0 

16 

27 

2 

1932.694 

1789.029 

55896.246 

0 

16 

27 

3 

1932.528 

1788.887 

55900.680 

0 

17 

27 

1 

1929.505 

1786.296 

55981.766 

0 

17 

27 

2 

1929.042 

1785.899 

55994.203 

0 

17 

29 

1 

1934.401 

1790.492 

55850.578 

0 

17 

29 

2 

1933.912 

1790.072 

55863.672 

0 

17 

29 

3 

1933.647 

1789.845 

55870.754 
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The  dipole  matrix-elements  which  are  necessary  for  the  calculation  of  the  transition 
probabilities  were  evaluated  using  an  ab-initio  method  in  which  the  required  wave  func¬ 
tions  for  both  the  X3£g‘  and  B3SU'  states  were  obtained  from  a  multi-reference,  first- 
order  configuration-interaction  calculation  with  a  large  Slater  basis  set  [27].  The  vibrational 
wave  functions  for  both  states  were  obtained  from  a  numerical  solution  of  the  one¬ 
dimensional  Schrodinger  equation  using  the  RKR  potentials  tabulated  in  Ref.  18.  The 
resulting  oscillator  strengths  for  the  v"(0)  to  v’(l-12)  transitions  agreed  with  low- 
temperature  measurements  [28]  to  within  5%.  The  transition  line-widths  were  assumed  to 
be  dominated  by  the  predissociation  rates  of  the  upper  states  which  are  given  by  Ref.20. 
The  collision  broadening  was  included  using  the  following  empirical  equations  derived  by 
Cann  et  al.[29]: 


to,  -  0.3P  (13) 

where  the  pressure  P  is  given  in  atmospheres.  For  T=300K  and  P=latm  Ai/p«8.43GHz. 
The  Doppler  broadening  was  assumed  to  have  a  Gaussian  line  shape  with  a  width  of  [30]: 


to,  -cv4ls 


(14) 


Where  yo=1.68‘1015  Hz  is  the  exciting  frequency  and  M  =  32,1.67’10‘27Kg  is  the 
molecular  mass  of  oxygen.  At  room  temperature  the  Doppler  broadening  is 
Au0  =  3.67GHz. 

The  total  line-width  of  the  transition  is  determined  by  the  accumulation  of  the  Dop¬ 
pler  broadening,  the  pressure  broadening  and  the  predissociation  line-width.  The  Doppler 
broadening  is  comparable  in  magnitude  to  the  predissociation  line-width.  The  final  line- 
width  however  is  dominated  by  the  pressure  broadening  mechanism.  Both  the  pressure 
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and  predissociation  mechanisms  have  a  Lorentzians  line-shape,  while  the  Doppler 
broadening  has  a  Gaussian  line-shape.  A  convolution  between  the  Lorentzian  line-shape 
and  the  Gaussian  shape  must  be  made  to  include  the  Doppler  broadening  in  the  calcula¬ 
tion.  This  however  requires  resort  to  either  polynomial  techniques  or  to  a  numerical  ap¬ 
proach  for  the  calculation  of  the  convolution  integral.  To  simplify  the  calculation  we  as¬ 
sumed  that  the  shape  of  the  Doppler  broadening  is  also  a  Lorentzian.  This  introduces  a 
relatively  small  error  but  has  the  advantage  of  a  much  simplified  calculation  and  a  broader 
range  of  analysis.  The  combined  line-shape  is  assumed  therefore  to  be  a  Lorentzian  with  a 
line-shape  function  gfy)  [21]: 


g(*0 


where  the  total  line  -width  Aut  is: 


(15) 


Aut  M  Avu  +  Avp  +  Avd 


(16) 


3.1.2  The  Fluorescence  Signal 

The  observed  fluorescence  signal  is  expected  to  be  proportional  to  the  initial 
population  of  the  ground-state,  N10.  Fundamentally,  this  population  is  proportional  to  the 
gas  number  density,  N0,  and  is  also  a  function  of  the  temperature  dependent  Boltzmann 
distribution  of  the  molecules  among  the  rotational  states.  The  variation  with  temperature 
of  the  population  density,  Na  of  an  individual  rotational  state  may  be  described  by  the  fol¬ 
lowing  equation  [31]: 


30 


(17) 


_Ni0-ShcBv 

j  KT 


where  S  is  the  Honl-London  factor  h  is  the  Planck's  constant,  Bv  is  the  rotational  constant, 
E  is  the  energy  of  the  rotational  state,  K=  138- 10*23  J/°K  is  Boltzmann’s  constant  and  N0  = 
0.2P/KT  is  the  density  of  oxygen  in  air  at  a  pressure  P  and  temperature  T. 

In  most  experiments  there  exists  a  temperature  at  which  the  population  density  of 
the  ground  state  reaches  a  maximum.  For  illustration,  at  constant  pressure  conditions,  the 
temperature  for  maximum  population  can  be  calculated  using  equation  (17)  and  is: 

Tmax  =  E/2K  (18) 

It  is  evident  from  eqn.17  that  the  LIF  and  FD  processes  depend  both  on  the  gas 
temperature  and  the  density  which  will  need  to  be  measured  independently  before  the 
velocity  can  be  deduced  from  the  FD  signal.  We  are  presently  developing  an  ArF  LIF 
technique  for  the  measurement  of  the  temperature  of  atmospheric  oxygen  and  a  Raman 
technique  for  the  measurement  of  the  density  of  air.  Both  techniques  will  be  used  simul¬ 
taneously  with  the  Doppler  velocimetry  technique  and  therefore  will  provide  the  additional 
necessary  information. 

When  an  intense  laser  is  used  for  this  measurements,  a  large  fraction  of  the  ground 
state  molecules  is  excited.  Owing  to  the  predissociation  effect,  only  few  molecules  relax 
back  to  their  original  state.  Therefore,  the  ground-state  population  during  the  LIF  process 
may  depend  also  on  the  laser  fluence.  While  replenishment  of  ground-state  molecules  can 
occur  by  several  mechanisms,  one  can  discard  the  effects  of  rotational  radiative  transitions 
which  are  negligible  because  electric-dipole  transitions  between  rotational  states  in  the  02 
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molecule  are  forbidden  and  the  allowed  magnetic-dipole  transitions  are  extremely  weak 
[32].  Thus,  only  collisional  relaxation  may  redistribute  the  population  among  rotational 
states.  However,  although  the  time  between  collisions  at  standard  temperature  and  pres¬ 
sure  conditions  is  0.14  ns,  it  will  be  comparable  with  the  15  ns  duration  of  the  laser  pulse  at 
low  densities.  Therefore,  at  the  low  gas  densities  which  may  occur  in  hypersonic  air  flows, 
the  absorbing  states  which  were  radiatively  depleted  will  not  be  equilibrated  and  the  strict 
linear  relation  between  the  total  fluorescence  and  the  population  density  may  be  violated. 

To  calculate  the  population  of  the  absorbing  state  and  the  subsequent  fluorescence 
signal,  including  the  effects  of  predissociation,  the  absorption  and  the  fluorescence 
processes  are  described  by  the  following  rate  equations: 


d(N,/gl)  _  N2/g2  N,  Ng_ 

*  '  ^  ,2(g,  g, 

d(N2/g2)  _*N2/g2  N2/g2  Nj  N2 

~ - ; - -  +  ”12'  —  - 


(19a) 

(19b) 


where  W12=Ia12/hv,  I  is  the  laser  intensity,  c12  is  the  absorption  cross-section,  hv  is 
the  photon  energy,  N,  and  N2  are  the  population  densities  of  the  lower  and  upper  states 
respectively,  gt  and  g2  are  their  degeneracies,  r21  is  the  spontaneous  radiative  lifetime  for 
the  transition  from  the  upper  state  to  the  lower  state,  rf  is  the  fluorescence  lifetime  of  the 
upper  state  including  transitions  to  all  the  accessible  levels  in  the  ground  electronic  state, 
and  rD  is  the  predissociation  time.  For  all  the  transitions  which  can  be  excited  by  the  ArF 
laser,  rQ  is  of  the  order  of  10'11  sec  while  r12  and  rf  are  approximately  10'5  sec  and 
W12  =  107  sec'1.  Because  the  predissociation  rate  exceeds  all  the  other  rates,  we  can  assume 
that  N2  <  <  N1  and,  therefore,  neglect  the  first  terms  in  eqns.(19a)  and  (19b).  The  rate 
equations  can  then  be  decoupled  and  solved  independently  to  give  the  result 
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N,»N10exp(-W12t) 


(20a) 


and 


glN  =  W12N10exp(-W12t)  W12N10  exp(-t/rD) 

g2  2  !A D-Wi2  Vr0  •  W12 

By  assuming  that  W12<  <  1  /rD,  Eq.(19b)  may  be  reduced  to 

^N2»r0W,2N10exp(-W,2t) 

The  fluorescence  signal  may  now  be  defined  as  the  total  number  of  fluorescence 
photons  which  are  emitted  during  a  laser  pulse  of  duration,  by  the  excited  molecules 
enclosed  in  a  volume  which  is  defined  by  the  beam  cross-section.  A,  and  the  length  of  the 
observed  beam,  L.  It  is  expressed  as 


(20b) 


(20c) 


S,  =  A-Lf  ?!  dt 

rf 

=  AL  JL  N,0(-|)  [l-exp(-W,2rp)]  (21) 

where  the  integration  has  been  simplified  by  assuming  that  the  spatial  and  temporal  inten¬ 
sity  distributions  of  the  laser  are  uniform.  With  those  assumptions,  Eq.  (21)  may  be  rewrit¬ 
ten  in  terms  of  the  laser  pulse  energy,  by  describing  the  average  intensity  as: 

1=  5e_  (22) 
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The  fluorescence  signal  is  then: 


s,  =  A'L(— rN0'(— )[l-exp(-(  Mil)]  (23) 

rf  1U  g,  A‘hi/ 

Note  that  this  expression  is  a  non-linear  function  of  the  laser  fluence,  Ep/A,  which,  in  turn, 
determines  the  extent  of  the  depletion  of  the  ground  state  population  by  the  predissocia¬ 
tion  process.  However,  when  W12rp  <  <  1,  or  when  the  rate  of  collisional  redistribution  is 
sufficiently  high,  the  depletion  of  the  ground-state  population  may  be  assumed  to  be  negli¬ 
gible.  The  fluorescence  signal  then  becomes  a  linear  function  of  the  laser  intensity  accord¬ 
ing  to: 


Sf(T»  =  L  A  •”»-  N,0(T,P)(ii)  (24) 

phv  rf  g1 

where  the  dependence  of  the  fluorescence  signal  on  the  temperature  and  density  of  the  gas 
volume  being  probed  is  indicated. 


3.13  UNCERTAINTY  IN  THE  VELOCITY  MEASUREMENT 

A  primary  goal  of  this  study  is  to  estimate  the  uncertainty  in  the  velocity  measure¬ 
ments.  For  this  analysis,  the  uncertainty  in  the  velocity  measurement  is  limited  by  the  error 
in  the  fluorescence  measurement  which  in  turn  is  assumed  to  be  limited  by  the  photon- 
statistical  noise  (the  shot-noise).  In  addition,  the  fluorescence  signal  for  the  temperature 
measurement  and  in  the  simultaneous  Raman  signal  used  for  the  required  measurement  of 
density  are  also  assumed  to  be  limited  by  the  shot-noise.  These  assumptions  were  verified 
by  us  experimentally  [2].  The  statistical  nature  of  any  photoemission  or  photodetection 
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process  results  in  fluctuations  in  the  total  number  of  photoelectrons,  N#,  which  are 
produced  in  the  photodetector  owing  to  the  fluorescence  from  each  laser  pulse.  These  fluc¬ 
tuations  appear  as  noise  in  the  fluorescence  signal  and  are  characterized  by  the  expression 
[33] 

ANe  =  VNe  (25) 

To  estimate  the  number  of  photoelectrons,  the  parameters  of  the  collection  optics 
must  be  included  in  the  calculation.  For  example,  in  a  typical  experiment  the  fluorescence 
from  a  sample  volume  of  length,  L  (taken  to  be  1  mm),  is  collected  by  a  system  of  collec¬ 
tion  optics  with  a  collection  efficiency,  r?c  (equal  to  .015  for  f/2  optics)  and  detected  by  a 
photomultiplier  with  a  quantum  efficiency,  tjd  (typically  equal  to  .15).  The  detection  sys¬ 
tem  may  also  contain  a  bandpass  filter  with  a  transmission,  T#,  which  is  used  to  reject  back¬ 
ground  radiation  emitted  by  fluorescing  optical  components.  The  result,  which  is  a  function 
of  the  gas  temperature.,  T,  and  density,  p,  is  then: 

N,frj>)“Sf(Tj>)»,ndT(L  (26) 

where  Sf  is  obtained  from  Eq.(24).  While  temperature  and  density  are  deduced  from  a 
direct  fluorescence  or  Raman  intensity  measurements,  the  velocity  depends  on  the  dif¬ 
ference  between  two  fluorescence  measurements,  N„.  The  relative  velocity  measurement 
uncertainty’  can  then  be  evaluated  from 
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where  A p  is  the  uncertainty  in  the  independent  density  measurement  which  is  based  on 
the  photon-statistical  noise  associated  with  the  Raman  signal  and  AT  is  the  uncertainty  in 
the  temperature  measurement.  The  sensitivity  of  the  fluorescer.ce  difference  signal  to 
velocity  is  determined  by  the  derivative  term  in  the  denominator  of  Eq.(27).  When  the 
photon-statistical  noise  is  high,  the  uncertainty  can  be  reduced  by  selecting  transitions  with 
high  sensitivity,  e.g.,  transitions  with  a  narrow  line-width  as  proposed  here. 

3.2.  RESULTS 

3.2.1  Selection  of  the  Optimal  Absorption  Line 

The  excitation  spectrum  of  oxygen  is  obtained  by  exciting  the  gas  with  a  narrow 
band  laser  -  which  is  tuned  over  the  absorption  bands  within  its  tuning  range-  while  observ¬ 
ing  the  broadband  fluorescence  signal.  The  excitation  spectrum  originating  from  the  ab¬ 
sorbing  transitions  from  v"(0)  to  v’(13),  v’(14),  '/(IS),  v’(16)  and  to  v’(17)  was  calculated  by 
our  model  while  assuming  that  the  ArF  laser  is  continuously  tunable  between  192.9  nm  and 
193.7  nm.  Consequently,  the  Raman  shifted  beam  in  the  calculation  was  tunable  between 
178.59  nm  and  179.27  nm.  The  fluorescence  signal  was  assumed  to  be  normalized  by  the 
laser  pulse  energy  which,  in  turn,  is  measui  ed  from  a  portion  of  the  laser  beam  going  to  the 
sample  volume.  The  beam  paths  through  room  air  to  the  sample  volume  and  to  the  pulse- 
energy  detector  had  identical  lengths.  Therefore,  the  normalization  removed  the  effects 
both  of  pulse-to-pulse  energy  variations  and  of  absorption  by  room  air  which  varies  with 
laser  wavelength. 

Figure  6  presents  the  predicted  normalized  excitation  spectrum  of  room  air  which  is 
excited  by  a  Raman  shifted  ArF  laser  with  a  line  width  of  9GHz  and  an  energy  of  1 
mJ/pulse.  The  calculated  total  fluorescence  was  multiplied  by  a  transmission  function 
•‘'hich  is  centered  at  240  nm,  has  a  bandwidth  of  10  nm  and  a  peak  transmission  of  20%. 
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This  transmission  function  simulates  the  filtration  which  is  necessary  for  the  removal  of  the 
radiation  which  is  scattered  at  the  the  laser  frequency  and  of  fluorescence  which  may  be 
induced  by  the  wind  tunnel  windows.  The  abscissa  is  marked  with  the  fundamental 
wavelengths  of  the  laser.  Most  of  the  lines  which  can  be  observed  in  this  spectrum  are  com¬ 
posed  of  overlapping  transitions.  Only  the  v’(0)  to  ^(13)  R13  transition  at  193.13  nm  and 
the  v"(0)  to  v’(13)  P15  transition  at  193.44  nm  appear  sufficiently  isolated  from  the  other 
transitions.  However,  even  these  transitions  are  composed  of  three  fine  structure  lines 
which  are  due  to  the  spin-orbit  splitting.  These  lines  are  characterized  by  the  total  angular 
momentum  J  where: 


1  =  1: 

J  =  K-1 

1=2: 

ii 

* 

1=3: 

J  =  K+ 1 

and  where  K  is  the  total  quantum  number  apart  from  spin.  Figure  7  presents  the  detailed 
spectrum  of  the  P15  line.  The  observed  structure  is  the  convolution  of  the  laser  line  with 
each  of  the  fine  structure  absorption  lines  which  are  labeled  by  their  I  values.  As  a  result 
of  this  convolution  the  1=2  line  which  is  enclosed  between  the  1  =  1  and  1=3  lines  appears 
stronger  than  the  two  other  lines. 

Figure  8  presents  the  variation  with  gas  temperature  of  the  total  fluorescence  which 
is  observed  following  the  selective  excitation  of  each  of  the  fine  structure  lines  in  the  P15 
transition.  The  maximum  fluorescence  is  observed  at  270  K.  This  is  also  the  temperature  at 
which  the  population  of  the  ground  state  peaks.  The  temperature  for  maximum  fluores¬ 
cence  for  1=2  which  is  predicted  by  equation  (18)  is  249  K.  The  slight  deviation  from  the 
the  prediction  of  equation  (18)  can  be  attributed  to  the  contribution  to  the  total  fluores¬ 
cence  of  other  transitions  which  coincide  with  the  wings  of  the  laser  spectral  distribution. 
An  error  of  Wo  in  the  prediction  of  the  total  fluorescence  can  account  for  an  error  of  10  K 
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Fig. 8.  Tlie  variation  with  air-temperature  of  the  total  fluorescence  which  is  observed  following  the  selec 
tive  excitation  of  each  of  the  fine  structure  lines  of  the  v"(0)-v*(13)  PI 5  transition. 


in  the  prediction  of  the  temperature  for  maximum  fluorescence.  When  the  temperature  is 
below  400  K,  fluorescence  emitted  by  the  transition  from  1=2  is  the  highest.  However,  as 
temperature  increases  the  intensity  of  the  vH(0)  to  ^(17)  P29  1  =  1  line  which  overlaps  the 
1=3  line  increases.  As  a  result,  the  intensity  of  the  fluorescence  which  is  excited  when  the 
laser  is  tuned  to  the  1=3  line  increases  faster  than  the  intensity  of  the  other  two  lines  in  the 
P15  manifold.  When  there  is  an  overlap  between  two  transitions  which  originate  from  two 
rotational  levels  with  different  N  quantum  numbers,  the  analysis  of  the  temperature  de¬ 
pendence  of  the  fluorescence  must  include  the  contribution  of  both  lines.  Therefore,  to 
simplify  the  analysis,  only  the  R13  line,  which  does  not  overlap  any  other  transition,  should 
be  used  for  these  measurements. 

Figure  9  presents  the  temperature  dependence  of  the  total  fluorescence  which  is 
emitted  when  the  laser  is  tuned  to  the  v"(0)  v’(13)  R13  transitions.  All  the  fine  structure 
lines  behave  similarly.  The  temperature  for  maximum  fluorescence  when  1  =  1  is  lower  than 
the  temperature  for  maximum  fluorescence  for  1=3.  This  is  predicted  by  eqn.(18).  It  is 
seen  from  this  figure  that  the  curves  which  represent  the  three  fine  structure  lines  behave 
similarly.  This  indicates  that  the  contribution  of  other  transitions  to  the  total  fluorescence 
of  each  of  the  fine  structure  lines  is  negligible  for  all  temperatures  up  to  1500  K.  With 
the  availability  of  a  single  absorbing  line,  the  measurement  technique  will  resemble  the  ap¬ 
proach  of  Hiller  et  al.[16]  and  Hassa  et  al.  [17].  The  intensity  of  the  fluorescence  induced 
by  the  laser  in  molecules  moving  opposite  to  the  beam  will  be  subtracted  from  the  fluores¬ 
cence  intensity  induced  in  molecules  moving  in  the  direction  of  the  beam  propagation.  In  a 
typical  Differential  Intensity  Measurement  experiment,  a  laser  pulse  will  be  guided  into  the 
wind-tunnel  at  a  slant  angle  with  the  axial  flow  direction.  This  allows  detection  of  the  Dop¬ 
pler  shift  associated  with  the  main  velocity  component.  After  crossing  the  wind-tunnel  and 
after  a  delay  which  exceeds  the  pulse  duration,  the  beam  will  be  retro-reflected.  The  delay 
is  necessary  for  the  separation  of  the  fluorescence  induced  by  the  incoming  beam  from  the 
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Fig.9.  The  variation  with  air-temperature  of  the  total  fluorescence  which  is  observed  following  the  selec 
tive  excitation  of  each  of  the  fine  structure  lines  of  the  v"(0)-v’(13)  R13  transitions. 


fluorescence  induced  by  the  retro-reflected  beam.  The  difference  in  intensity  between  the 
two  fluorescence  bursts  can  be  related  theoretically  or  experimentally  [17]  to  the  flow 
velocity.  Figure  10  presents  schematically  the  interaction  of  the  laser  beam  with  the  R13 
line.  The  iaser  is  tuned  to  the  short  wavelength  (high  frequency)  side  of  the  1=3  line.  The 
excitation  spectra  which  are  expected  for  the  Doppler  up-shifted  and  Doppler  down-shifted 
frequencies  were  marked  by  the  dashed  lines.  The  intensity  differences  between  t.?  i 
fluoi  escence  induced  by  the  Doppler  up-shifted,  the  Doppler  down-shifted  and  the  un¬ 
shifted  frequencies  are  evident  in  this  figure. 

Figure  11  presents  the  photon  count  difference,  ND,  i.e.  the  difference  between  the 
fluorescence  intensity  following  the  excitation  of  the  the  v"(0)  to  v’(13)  R13  transition  by  a 
laser  beam  which  is  propagating  opposite  to  the  flow  and  between  the  fluorescence  which  is 
excited  by  a  laser  which  co-propagates  with  the  flow,  vs.  the  laser  wavelength  while  the 
velocity  is  selected  as  a  free  parameter.  The  figure  represents  velocities  ranging  from  100 
m/sec  to  2000  m/sec  in  intervals  of  100  m/sec.  The  upper  velocity  limit  is  approximately 
the  highest  air  speed  which  can  be  obtained  for  a  stagnation  temperature  of  1500  K..  For 
each  of  the  fine  structure  lines,  the  difference  is  negative  when  the  laser  is  tuned  to  the 
high  frequency  side  of  the  line  and  is  positive  when  the  laser  is  tuned  to  the  low  frequency 
side.  This  figure  can  be  used  to  select  the  laser  wavelength  which  will  induce  the  FD  which 
is  the  most  sensitive  to  air  velocity.The  maximum  difference  is  expected  at  the  inflection 
point  of  the  line.  Therefore,  two  maxima  are  expected  for  each  of  the  fine  structure  lines. 
The  measurement  of  the  difference  in  the  fluorescence  is  most  sensitive  to  the  flow  velocity 
when  the  laser  is  tuned  to  193.1272  nm  which  is  at  the  high  frequency  side  of  the  1=3  line 
and  when  it  is  tuned  to  193.1392  nm  which  is  on  the  low  frequency  slope  of  the  1  =  1  line. 
Tuning  the  laser  to  the  center  of  a  single  absorption  line  will  result  in  no  velocity  sensitivity. 
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3.2.2.  The  Effect  of  Temperature  and  Density  on  the  Measurement  Accuracy 

The  measurement  of  N0  is  expected  to  vary  with  temperature  and  density  very  much 
like  the  variation  of  the  total  fluorescence  (eqn.  17).  Therefore,  the  measurement  of  the 
velocity  using  this  technique  must  rely  on  the  independent  measurement  of  the  tempera¬ 
ture  and  the  density. 

The  dependence  of  the  FD  measurement  on  temperature  and  density  may  be 
reduced  by  dividing  the  FD  with  the  total  fluorescence,  Nc,  induced  in  the  absence  of  a 
Doppler  shift.  This  signal  contains  the  same  temperature  and  density  information  as  the 
FD  but  is  independent  of  the  flow  velocity.  This  additional  fluorescence  measurement  can 
be  obtained  by  using  a  third  laser  beam  which  propagates  perpendicularly  to  the  flow.  In 
cases  where  the  slope  of  the  wing  of  an  absorption  line  can  be  approximated  as  a  straight 
line,  this  correction  of  the  FD  will  eliminate  the  temperature  and  density  dependence.  In 
most  cases  however,  such  an  assumption  may  not  be  warranted  and  therefore,  this  correc¬ 
tion  is  expected  only  to  reduce  the  temperature  and  density  effects  on  the  FD  signal. 

The  variation  of  N„  with  velocity  was  calculated  while  selecting  the  temperature  as  a 
free  parameter.  The  laser  wavelength  was  X  =  193.1272  nm.  This  the  the  wavelength  at  the 
high  frequency  side  of  the  R13/I3  line  at  which  the  velocity  sensitivity  is  the  highest.  Figure 
12  presents  the  results  of  this  calculation  for  temperatures  ranging  from  100  K  to  1500  K. 
This  is  the  range  of  interest  in  most  hypersonic  wind  tunnel  applications.  The  count  dif¬ 
ference  Nd  in  this  figure  increases  with  temperature  when  the  gas  temperature  is  below  250 
K  and  than  decreases  as  temperature  increases.  This  is  similar  to  the  variation  of  the  total 
fluorescence  with  temperature  as  predicted  by  figure  9.  In  addition,  the  slope  of  all  curves 
decreases  at  high  velocities.  Therefore,  the  measurement  sensitivity  is  expected  to  decrease 
as  the  air  velocity  approaches  2000  m/sec.  It  is  apparent  from  this  figure  that  the  measure¬ 
ment  of  N0  can  not  be  correlated  with  velocity  without  an  independent  measurement  of  the 
temperature.  Furthermore,  according  to  eqn.(17),  the  population  of  the  ground  state  is 
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proportional  to  the  total  density  and  the  absorption  cross-section  varies  with  pressure  due 
to  the  pressure  broadening,  Therefore,  Nc  aiso  varies  with  density  or  pressure  and  the 
velocity  measurement  must  also  include  a  separate  density  measurement.  As  a  result,  the 
uncertainty  in  the  measurement  of  the  velocity  will  be  affected  by  the  uncertainties  in  the 
measurements  of  the  density  and  the  temperature. 

To  evaluate  the  effectiveness  of  the  temperature  and  density  correction  technique, 
we  repeated  the  calculation  of  the  FD  and  used  the  total  fluorescence  induced  by  the  un¬ 
shifted  beam  for  normalization.  Figure  (13)  presents  the  variation  with  velocity  of  the  tem¬ 
perature  and  density  corrected  FD,  ND = N0/Nc.  Several  temperatures  were  selected  as 
free  parameters  while  the  pressure  was  kept  constant  at  1  atmosphere.  The  other  condi¬ 
tions  for  this  calculation  are  identical  to  those  described  for  figure  (12).  This  calculation 
presents  the  temperature  dependence  of  the  corrected  FD.  It  is  seen  that  for  low  tempera¬ 
ture,  when  the  contribution  of  both  the  pressure  and  temperature  broadening  to  the  total 
line  width  varies  rapidly  the  corrected  FD  is  sensitive  to  temperature.  However,  above 
500K  the  corrected  FD  is  nearly  independent  of  temperature.  This  strong  temperature  de¬ 
pendence  below  500  K  was  attributed  to  the  rapid  changes  in  the  pressure  and  Doppler 
broadening  (eqns.  (13)  and  (14)).  While  the  pressure  line-width  decreases  with  the  tem¬ 
perature  the  Doppler  line-width  increases.  Above  500  K,  The  sum  of  these  two  mechanisms 
was  almost  independent  of  temperature.  This  is  where  the  temperature  correction  becomes 
effective. 

To  evaluate  the  density  dependence  of  the  corrected  FD  we  repeated  the  calcula¬ 
tion  with  the  pressure  as  the  free  parameter.  Figure  14  presents  the  variation  with  velocity 
of  the  temperature  and  density  corrected  FD.  Here  the  temperature  was  kept  at  100K.  The 
sensitivity  of  the  FD  to  pressure  was  found  to  be  the  greatest  at  this  temperature.  It  is  seen 
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that  although  the  FD  has  been  corrected  by  the  normalization  it  still  presents  a  strong  den¬ 
sity  dependence.  The  error  term  (5N0/ap)VJAp  which  is  introduced  into  the  velocity 
measurement  due  to  uncertainty  in  the  density  measurement  was  found  to  increase  with 
density.  Again,  the  increasing  pressure  line-width  was  found  to  be  the  cause  for  the  density 
dependence  of  the  corrected  FD. 

In  a  previous  work  [2],  we  determined  the  uncertainty  which  is  associated  with  the 
measurement  of  the  temperature  and  density  in  hypersonic  wind  tunnel  air  flows  using  an 
LIF  and  Raman  technique  and  assuming  that  the  measurement  will  be  limited  by  the 
photomstatistical  noise.  The  uncertainty  in  the  temperature  measurement,  AT/T,  when 
the  air  density  is  0.01  amagat  and  the  laser  energy  at  193.4  nm  is  15  mJ/pulse  was  found  to 
be  less  than  0.5%.  The  uncertainty  in  the  density  measurement,  A  p/p,  using  a  15  mJ  pulse 
to  excite  Raman  scattering  from  0.01  amagat  of  air  was  found  to  be  less  than  2.7%.  This  is 
however  a  conservative  estimate.  The  pulse  energy  at  193.4  nm  in  most  applications  may 
exceed  100  mJ.  Therefore,  the  uncertainty  in  the  density  measurement  is  expected  to  be 
less  than  1%.  The  model  may  now  be  applied  to  estimate  the  uncertainty  v/ith  which 
velocity  can  be  measured  using  an  achievable  laser  performance  at  the  gasdynamic  condi¬ 
tions  of  interest.  Examples  calculated  using  Eq.  (27)  are  based  on  the  assumption  that 
photon-statistical  noise  is  the  limiting  cause  for  the  uncertainty  in  the  velocity  measure¬ 
ment.  However,  we  shall  a*so  show  examples  in  which  the  minimum  Noise  to  Signal  Ratio 
(NSR)  is  limited  to  1%  to  simulate  the  case  where  the  noise  is  induced  by  Electromagnetic 
effects  or  by  digitization.  The  uncertainty  of  an  LIF  temperature  measurement  was  calcu¬ 
lated  using  Eq.(27)  for  gas  pressures  of  .01  atmospheres  and  temperatures  varying  from  100 
K  to  1500  K.  This  covers  the  range  of  most  applications.  For  example,  at  a  free-stream 
Mach  number  of  10,  a  total  pressure  of  100  atm.,  and  a  total  temperature  of  1000  K,  the 
flow  field  over  a  typical  hypersonic  aerodynamic  body  has  densities  which  vary  from  0.01  to 
0.05  amagat  and  temperatures  from  50  to  250  K.  The  stagnation  temperature  however,  may 
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approach  1000  K.  The  calculations  were  made  assuming  that  the  Raman  shifted  laser  pulse 
energy  is  1  mJ,  with  fluorescence  and  Raman  scattering  observed  from  a  1-mm  beam  path 
with  i/2  collection  optics  and  a  detector  with  15%  quantum  efficiency.  These  optical  condi¬ 
tions  are  equivalent  to  a  laser  fluence  of  .1  J/cm2  in  a  sample  volume  with  a  beam  cross  ' 
section  of  1  mm2.  The  uncertainty  in  the  temperature  measurement  was  assumed  to  be 
.05%  and  the  limiting  uncertainty  in  the  density  measurement  was  assumed  to  be  1%.  The 
density  error  term  (3N0/5p)v  TAp  was  evaluated  using  the  calculation  of  fig.  (14)  for 
P  =  1  atm  and  T  =  100  K.  This  is  the  largest  error  in  tills  range.  The  temperature  error  term 
was  calculated  separately  for  each  temperature. 

The  results  are  shown  in  figure  (15)  for  the  laser  tuned  to  the  high  frequency  slope 
of  the  13/R13  rotational  transition  in  the  v"(0)  -  v’(4)  vibrational  band.  The  solid  line  is 
based  on  the  assumptions  that  the  uncertainty  in  the  velocity  measurement  is  composed  of 
the  uncertainty  in  the  temperature  and  density  measurements.  However,  unlike  the  as¬ 
sumptions  in  eqn.27  the  uncertainty  in  the  FD  was  assumed  to  be  limited  to  1%.  Therefore, 
the  solid  line  presents  the  most  conservative  estimate  of  the  measurement  uncertainty.  It 
indicates  however  that  at  the  most  adverse  conditions  the  measurement  uncertainty  may  be 
kept  below  10%.  The  next  estimate  presents  the  uncertainty  associated  with  a  velocity 
measurement  using  the  temperature  and  density  corrections  technique.  Here,  the  uncer¬ 
tainty  can  be  kept  below  5%.  The  dotted  line  shows  the  results  if  the  minimum  NSR  in  the 
FD  is  limited  to  1%  by  EMI  or  other  noise  sources  but  the  temperature  and  density  are 
determined  precisely.  This  presents  the  uncertainty  in  the  velocity  measurement  which  may 
be  obtained  in  some  experimental  conditions.  Finally,  the  dashed  line  presents  the 
measurement  uncertainty  when  the  FD  is  limited  by  the  shot  noise  and  the  temperature 
and  density  are  determined  precisely.  This  curve  presents  the  physical  limit  for  the  velocity 
measurement  uncertainty  using  this  Doppler  velocimetry  technique.  It  is  seen  that  the  ac¬ 
curacy  of  the  velocity  measurement  using  Oz  as  the  tracer  molecule  can  not  exceed  0.1% 
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Fig.  15.  Velocity  measurement  uncertainty  Av/v  vs.  temperature  for  the  laser  tuned  to  the 
v’(0)-v’(13)  R13  transition.  The  solid  line  presents  the  total  uncertainty  which  is 
composed  of  the  uncertainty  in  the  temperature  and  density  measurements  and 
aNd//N0  =  1%.  The  presents  the  uncertainty  when  using  the  temperature  and 
density  corrections.  The  dotted  line  presents  the  uncertainty  when  the  only  error  is 
AN0/N0  =  1%.  The  dashed  line  presents  the  measurement  uncertainty  when  the 
only  error  is  N0  =  v/N0. 


4.  CONCLUSIONS 


A  new  velocimetry  technique  which  uses  a  single  pulse  from  an  ArF  laser  for  the 
generation  and  the  subsequent  detection  of  CO  tracers  was  analyzed  theoretically  and  ex¬ 
perimentally.  The  intensity  emitted  by  fluorescing  carbon  atoms,  which  are  the  dissociation 
products  of  atmospheric  C02,  was  expected  to  be  sensitive  to  air  velocity. 

An  alternative  technique  which  uses  the  Doppler  shift  experienced  by  moving  02 
molecules  for  the  measurement  of  their  velocity  was  evaluated  theoretically.  The  Doppler 
shift  is  measured  by  using  a  differential  intensity  measurement  technique.  In  this  approach 
two  counter-propagating  laser  beams  are  tuned  to  the  high  frequency  wing  of  one  of  the 
absorption  lines  of  02  molecules  within  the  air  flow.  In  the  absence  of  flow,  both  beams 
encounter  the  same  absorption  cross-section  and  induce  the  same  fluorescence  intensity. 
However,  when  the  absorbing  molecules  are  moving  along  the  laser  beams  together  with 
the  gas  flow,  the  absorbed  radiation  is  Doppler  shifted:  the  frequency  of  the  radiation  of 
the  beam  which  propagates  opposite  to  the  flow  is  Doppler  up-shifted  and  is  therefore  fur¬ 
ther  away  from  the  absorption  line-center.  Therefore,  fluorescence  intensity  induced  in 
molecules  which  move  opposite  to  the  beam  decreases.  The  other  laser  beam  -  which 
propagates  in  the  direction  of  the  N2  flow  -  is  experienced  by  the  moving  molecules  as  Dop¬ 
pler  down-shifted.  Therefore,  the  fluorescence  induced  by  this  beam  will  be  more  intense. 
The  difference  between  the  intensities  of  these  two  LIF  processes  is  sensitive  to  the  air 
speed. 

The  results  of  our  study  yielded  the  following  conclusions: 

When  the  air  temperature  is  below  2000  K  and  the  pressure  is  between  0.01  atmos¬ 
phere  and  10  atmospheres  the  fluorescence  emitted  by  the  dissociation  products  of 
C02  is  independent  of  temperature.  Above  this  temperature  the  equilibrium  density 
of  COz  decreases  rapidly. 
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The  total  fluorescence  emitted  by  the  dissociation  products  of  COz  strongly  depends 
on  the  energy  distribution  across  the  laser  beam.  Fluctuation  in  the  energy  distribu¬ 
tion  may  introduce  an  uncertainty  in  the  velocity  measurement  of  up  to  40%. 

A  preliminary  experiment  tested  the  sensitivity  of  the  total  fluorescence  which  is  in¬ 
duced  by  an  ArF  laser  in  the  dissociation  products  of  COz  to  supersonic  flows.  No 
variation  in  the  total  fluorescence  could  be  detected  for  air  speeds  up  to  Mach  2.2. 
It  is  believed  that  multiphoton  ionization  and  quenching  processes  interrupt  the 
photodissociation  of  C02  and  CO  and  the  UF  process. 

The  v"(0)  -  v’(13)  R13  transition  is  the  only  line  which  is  resonant  with  the  Raman 
shifted  ArF  laser  and  is  also  isolated  from  other  absorption  lines  in  the  Schumann- 
Runge  band. 

It  was  shown  theoretically  that  even  if  the  FD  signal  is  normalized  by  the  total 
fluorescence  which  is  emitted  by  stationary  molecules  the  normalized  signal  still 
presents  a  strong  density  and  temperature  dependence. 

When  the  FD  signal  is  not  corrected  for  density  and  temperature,  the  uncertainty  in 
the  velocity  measurement  for  all  temperatures  below  1500  K  and  for  densities  below 
1  atmosphere  is  less  than  10%. 

When  the  FD  signal  is  normalized  by  the  total  fluorescence  emitted  by  stationary 
molecules  the  uncertainty  in  the  velocity  measurement  for  all  temperatures  below 
1500  K  and  densities  below  1  atmosphere  is  less  than  5%. 

The  ultimate  accuracy  of  this  Doppler  velocimetry  technique  is  obtained  when  the 
temperature  and  density  are  accurately  known  and  when  the  fluorescence  is 
measured  within  an  uncertainty  which  is  limited  by  the  shot-noise.  The  uncertainty 
in  the  velocity  measurement  is  than  less  than  0.1%. 


55 


5. 


RECOMMENDATIONS 


1.  The  single  pulse  carbon  velocimetry  did  not  appear  to  have  the  necessary  sensitivity 
to  air  speed.  Therefore,  we  do  not  recommend  its  application  for  wind  tunnel 
diagnostics. 

2.  We  recommend  the  application  of  the  Doppler  velocimetry  in  oxygen  for  wind  tun¬ 
nel  diagnostics. 

3.  We  recommend  to  use  the  Fluorescence  Difference  (FD)  approach  for  single  pulse 
velocity  measurements. 

4.  We  recommend  to  use  the  absorption  in  the  v"(0)-v’(13)  R13  transition  for  these 
velocity  measurements. 

5.  We  recommend  to  normalize  the  FD  signal  by  the  total  fluorescence  induced  in  sta¬ 
tionary  molecules.  This  normalization  will  reduce  the  temperature  and  density  de¬ 
pendence  of  the  FD  measurement. 
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